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Multi Higgs doublet models are interesting extensions of the Standard Model that 
can be related to flavor. The reason is that most flavor models usually involve the 
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presence of several additional scalar fields. In this work we present an analysis that 
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shows that for renormalizable flavor models based on the cyclic group of order N, 
if there is one flavored SU(2) double Higgs per generation, the smallest N that can 
be used to reproduce the Nearest-Neighbor-Interaction texture for the quark mass 
matrices is N = 5. Results for the Higgs spectrum and consistency under K — K 



mixing in a specific model with are also presented. 



The Yukawa matrices of the Standard Model (SM) parametrize our ignorance on the pos- 
sible relations and perhaps origin of the mass spectrum and mixing angles of all fundamental 
fermions. Given the experimental determination of fermion masses (so far only mass squared 
differences in the case of neutrinos), the entries of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix (including its CP violating phase), and the values of the neutrino mixing angles en- 
coded into the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, it is possible to determine, 
in a model independent way, certain textures for the fermion mass matrices (obtained from 
the Yukawa matrices) that upon diagonalization reproduce those results. In particular, con- 
sidering for the moment the quark sector, the so-called Nearest-Neighbor-Interaction (NNI) 



texture 



2[ successfully reproduces the quark masses and CKM angles and phase. The NNI 
texture has the general form 

(1) 
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and can always be obtained from general mass matrices through a proper choice of flavor basis 
in the SM [ll,!^]- The lepton sector can be included by either assuming a diagonal mass matrix 
for the charged leptons and diagonalizing the neutrino mass matrix with the PMNS matrix, or 
by considering the possibility that the PMNS matrix receives contributions from both sectors. 
Extending the SM in order to incorporate descriptions of fermion masses and mixing angles 
thus necessarily requires either to predict or accommodate those specific textures. 

Intimately related to these issues is the scalar sector of the SM and its extensions. The 
spontaneous breaking of Electroweak (EW) symmetry in the SM, through the non-zero vacuum 
expectation value (vev) of a Higgs field, is what generates the mass matrices from the Yukawa 
matrices. The SM incorporates this mechanism in a minimal fashion, i.e. by postulating 
the existence of one doublet (under SU(2)) scalar field. Other choices and possibilities have 
been explored and there is a vast literature associated to them. Among the most popular 
extensions are those that consider the case (within many different settings) with two Higgs 
doublets, generically called two Higgs doublet models (THDM - see {^J for a recent review). 
Higgs triplets (in combination with doublets) have also been studied extensively In most 

cases these extensions are motivated by the scalar sector phenomenology and do not necessarily 
add much to the problem described above. Most cases involve flavor-blind Higgs fields that do 
not directly participate into the structure of the mass matrices. 

An interesting scenario where the Higgs fields necessarily participate in the flavor structure 
is the one of renormalizable flavor models (3, sj|. In this type of models a horizontal flavor 
symmetry, continuous or discrete, is added to the SM gauge group symmetry in such a way 
as to reproduce the observed mass and mixing angle patterns by only using renormalizable 
Lagrangians. This requirement has two immediate and interesting consequences: i) there must 
be more than one SU(2) doublet scalar and ii) at least some of them must transform non-trivially 
under the flavor symmetry. 

In this letter we consider the case of renormalizable models involving three Higgs SU(2) 
doublets (Hi, H2, H3), the SM gauge and fermion particle content, and a discrete Abelian 
flavor symmetry Z^. The purpose is to find the smallest that can be used to obtain the 
NNI textures in Eq.((T]) in a setting with three Higgses (the smallest realization of the NNI 
structure with two Higgs doublets requires Z4 qJ). This setting is interesting from the point 
of view of minimal extensions to the SM that could give useful hints to more elaborated and 
ambitious extensions. Studying the case of three SU(2) doublets is motivated by the fact that 
most flavor models, renormalizable and non-renormalizable, usually require the presence of 
several additional scalar fields. Our purpose here is to show that having one flavored Higgs per 
generation, can lead to interesting possibilities. 
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In order accomplish our goal, the following considerations are taken into account: 

i) Each non-zero entry in the mass matrices contains contributions from a single vev. The 
condition is required by the fact that if several Higgses contribute to the same entry, they 
necessarily will have the same charge, making them virtually indistinguishable. 

ii) One of the Higgses will be assigned neutral charge while the other two will be related 
by conjugation. This is necessary in order to use the same three Higgses in both the up 
and down quark sectors (in the case of even N it is also possible to consider one of the 
Higgses to have charge N/2 instead of the neutral charge, however this requires larger 
groups than the neutral case. See Appendix lAl). 

iii) One of the Higgses contributes exclusively to the 3 — 3 entry of the mass matrix for the 
up-type quark. This is motivated by the fact that the top quark is the heaviest and 
resembles the familiar case of models based non-Abelian flavor groups where the fermions 



are put in 2 © 1 representations 



10|. 



Given these considerations, the possible textures for the up-type quarks have the general 
form (k 7^ I ^ m) 
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where v^i >m ), (k, l,m) = 1, 2, 3 denote the Higgs vevs. 

The charge assignments for the fermion and scalar fields is parametrized as 



(2) 



Q ~ (91,92,93), U R ~ (ui,« 2) «3), Dr - (di,d 2 ,d 3 ) 
U = (#1, H 2 , H3) — (hi, h 2 , h 3 ) , 



(3) 
(4) 



where u { , di, hi G Z N , Q is the left-handed quark SU(2) doublet, Ur (D r ) is the up-type 
(down-type) right handed quark SU(2) singlet, and Hi = ia 2 H*. The charges of the 
bilinears in the Yukawa terms of the Lagrangian can be represented by 



yj- = qi + Uj mod(N). 



(5) 



We now present the results showing the complete analysis for one case. 
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Case Al: Without loss of generality, let k = 1, I — 2, and m = 3 for M\ x in Eq. 
From condition ii) above we have the following three possibilities for the Higgses charges: a) 
(hi = 0, h 2 = a, h 3 = —a), b) (hi = a,h 2 = 0, h 3 = —a), and c) (hi = —a, h 2 = a,h 3 = 0) with 
a G Z N . 

Case a) leads to the following constraints on the fermion bilinear charges Eq. 



y u u = 0, 3& = 0, 3& = -a, ^2 = 3% = a, (6) 



and 



^i u i, K 3 , y^(0.-M). (7) 

Using Eq. ([H]) we find that the fermion assignments become 

qi = -c, q 2 = -c- 3a, q 3 = -c - a, (8) 
u\ = c + 3a, u 2 = c, u 3 = c + 2a, (9) 
di = c + 3a, d 2 — c, d 3 = c + 2a, (10) 

where c, a G -Zat. The last step is to satisfy the relations in Eq. (j^) which for this case become: 



(3a, 2a, —3a) ^ < 





-a mod(N) . (11) 
a 



Since 2a = in Z 2 and 3a = in Z 3 , these two groups are discarded. For Z4, a G (0, 1, 2, 3). 
Since a 7^ 0, then a G (1, 2, 3). Since 2a 7^ 0, then a G (1, 3). If a = 1, then 3a = 3 = — 1 = —a 
and thus a ^ 1. must be discarded. For Z 5 a G (0, 1,2,3,4) and again, since a 7^ 0, then 
a G (1,2,3,4). If a = 1, then —a = 4, 2a = 2, 3a = 3, —3a = 2, and thus all conditions in 
Eq. fllip are satisfied and the smallest group that works in this case is Z 5 . 

A similar analysis for all other possibilities shows that 

• Z§ is the smallest group that can be used in this setting and corresponds to the case 
Al-a). Larger groups can also be used in this case. 

• Z 6 is the smallest possibility for case Al-c). Larger groups can also be used in this case. 

• All other cases do not satisfy the necessary conditions for any Z^. 

Thus, the smallest group that can be used within this setting is Z 5 . Note that the analysis 
includes the down-type quark sector as well by choosing the same charge assignments as those of 
the up-type. However, it is also possible to obtain the required texture with other assignments. 
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This is of course a model dependent issue. As for leptons, it is straightforward to generalize 
to the charged leptons by giving them the same charge assignments used for the down-type 
quarks. Neutrino masses and mixing can then be included either considering radiative mass 



generation 



11 



1311. n on-renormalizable terms, or seesaw through the introduction of right- 



handed neutrinos 14j . This is a model dependent question and we do not explore it further in 
this letter except to mention that for instance, in a model like the one in case Al-a), the lepton 

n 

sector can reproduce exactly the one presented in |8[ . A specific model and its phenomenology 
will be presented in a future publication. 

We note that, given the condition on the scalar charge assignments, it is possible to write 
a general invariant potential (a general discussion on Abelian symmetries in multi-Higgs 
models can be found in [ij]]). Since one of the Higsses is required to be neutral and the 
remaining two to be conjugate of each other, and for clarity, let us use the following notation: 
denote by H the neutral one and by $ OJ a = 1,2 the remaining two. The potential is then 
given by 

V(H, $ a ) = fi 2 \H\ 2 + + /4 (& a H + h.c) + /i? 2 ($I$ 2 + h.c) + A (\H\ 2 ) 2 

+ A Q (|$ a | 2 ) 2 + \ 0a \H\ 2 \<S> a \ 2 + A 12 |$i| 2 |$ 2 | 2 + ~X ab \^ b \ 2 + \' 0a & a HH^ a 
+ A 3 ($\H$lH + h.c\ , (12) 

where the terms proportional to /ioa and fi\2 are Zn soft breaking terms required in order to 
obtain the correct electromagnetic invariant vacuum. Denoting the scalars by 

, H+ I ' = [ , ^ 

-L(v + h + iAo) J \ (v a + <f) a + iAa) 

where Vq and v a denote the vevs of H and $ a respectively, the minimization conditions become 

2 Aq^^q + 2\ Z V 1 V 2 Vo + Ap^fo + Xoivlvp + A 2^0 + 2\ vl + 2V 1 /J,l 1 + 2v 2^02 nA \ 

N - 2v ' 1 j 

Mi = ^ , (15) 

2 A 2^ 2 + \3vrv + Ai 2 t> 2 t>2 + 2X2V2 + A 2t>2^ 2 + 2v 1 fi 2 2 + 2vqh1 2 

In order to perform a numerical analysis for the Higgs mass spectrum it is necessary to study 
a specific model. Taking as an example the case Al-a) for Z5 we have performed a scan over 
the parameter space of the model. We find that there are large regions of parameter space 
consistent with current experimental values and bounds. As an example we present in Table [I] 
one particular set of parameters (a complete study including a statistical analysis involving a 
X 2 fit, possible collider signatures, as well as a compelte analysis of the lepton sector, is under 
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preparation and will be published in a future paper) that gives the following spectrum: For the 
three CP-even scalar masses we obtain (in GeV): 125.7, 700.9, and 892.1; for the two massive 
CP-odd scalars we obtain (in GeV): 670.4 and 894.1; and for the two charged scalars (again in 
GeV): 678.7 and 895.3. 



vi (GeV) 


v 2 (GeV) 


^3 (GeV) 


/4 (GeV) 2 


& (GeV) 2 
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Af)2 
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0.63927 


-0.561199 


0.160189 


0.0779788 


-0.758485 


0.426743 


-0.543321 


-0.582515 


-0.0203623 



Table I: One set of parameter values consistent with current experimental data. 



Furthermore, since the model contributes to flavor changing transitions at tree level, we also 
incorporate the constraint coming from K — K mixing. This is done by computing AM^ from 



the effective Hamiltonian 



lfHial. 



«r=^EG W e W (it) 

i 

where GV = 1.16639 x 10~ 5 GeV -2 is the Fermi constant and C{ are the Wilson coefficients. 
Consequently, the K — K mixing is governed by the neutral scalar interactions with the first 
and second down type quark families, then 

Q™ = (sP L d)(sP R d) (18) 
Qf LL = (sP L d)(sP L d) (19) 
Qf™ = (sP R d)(sP R d) (20) 



3 

n LR _ I '"-dill's \ ST^ U2aUla , Q1 v 

where all fields are in the mass basis and U a b denotes the scalar mixing matrix. The K — K 
mixing is given by the off-diagonal term in the neutral .ff -meson mass matrix 

K _ AMk = x 10 _ 15 _ (24) 

12 M K y 1 



are 




16tt 2 / 
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where AM^ is given by 18 ] 



AM K = 2Re(K°\n^ s f 2 \K ) 



12tt 2 



M k F 2 kV2 B k [P 2 lr C^ r + Pf LL (Cf iL + Cf**)] , (25) 



where .FV = 160 MeV is the K-meson decay constant, = 497.6 MeV is the 



mass, and Bk and 772 include the QCD running effects. We follow the notation in 
pSLL = _ g _ 3) pLR = 30>6) V2 = a57; fi K = Q 85 ± o 15 Consequently, we obtain 

4 1 3 

^£ = -F 2 V2 B K (m d m s ) - V 

3 * m 



^-meson 
so that 



a=l 



yLR^aUia fySLL ( & 



'2a 



u 2 



(26) 



Taking into account that ~ 5 MeV and m s ~ 100 MeV, it is possible to establish that the 
constraint coming from KK mixing is satisfied in our model when 



J r c< 



3M* 



8.82 x 10- 10 (GeV 



AF^r]2B K m d m s 



(27) 



Taking the values in Table U we obtain 



J=l" = 5.5025 x lO^GeV" 4 . 



(28) 



In conclusion we have shown that the smallest cyclic group that can be used to generate the 
NNI textures for the quark mass matrices, in the context of three SU(2) Higgs doublets where 
one couples only to the 3 — 3 entry of the up-type quark mass matrix, is Z 5 . We have outlined 
the analysis that led to this result and have presented the general scalar potential for such a 
scenario. The results presented here can be of use to model builders interested in flavor models 
with three Higgs doublets. 
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Appendix A: Analysis for even N 



Recall from Eq. (j2J) all possible textures for the up-type quarks: 



' v k ^ 



Ml, ~ 



V k 

uj v rn 



1 y fc \ 



M U A2 ~ 



vi w fc 
vi v m 



' v k ^ 



w z vi 
Wfc w m 



' v k ^ 



v k w fc 
v t v rn 



1 y fc ^ 



M U B2 ~ 



vi Vi 
uj w m 



(Al) 



where v (k,i,m)i ^w) = 1, 2, 3 denote the Higgs vevs. 

In general the Higgses have charges H — {hi — a, /12 = b, h 3 
lets consider case A\ 

Ai) The fermion fields have the following charges: 

qi = a, q 2 = —c + a — a + 2b, q 3 = a — a + b 

Ui = 2a + c — a — 2b, u 2 = a — a, u 3 = a + 3b + c — a 

and the constraints become 



c) where a,b,c G Z^. Now 



(A2) 



(3^11,^22,^13,^31)^ < 



b mod(N) . 
c 



(A3) 



where, 



y u =2a + c- 2b, y 22 = -c + 2b, y 13 = a + 3b + c and ^31 = a - b + c. (A4) 

Taking into consideration the condition that two of the Higgses are related by conjugation, 
we have the following possibilities, a) H = (d,e,—e), b) H = (e,d,—e) and c) H = 
(— e, e, d) where, e,d G Z^. From the additional condition that the remaining Higgs is 
either neutral or, for N even, N/2, we have the two cases d = for any N and d = N/2 
for iV even. The case d = has been presented in the paper and we now present the 
analysis for d = N/2. In this case the constraints for each possibility of A\ are: 



a) 



y u =N- 3e, y 22 = 3e, y 13 = N/2 + 2e and y 31 = N/2 - 2e, 



(A5) 



which are satisfied with an Abelian symmetry group of order N > 8. 
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b) 

y u = e - N, y 22 = e + N, y 13 = 3N/2 and y 31 = -N/2, (A6) 
which can not be satisfied with any Abelian symmetry. 

c) 

y u = -4e + N/2, y 22 = -N/2 + 2e, y 13 = N/2 + 2e and y 31 = N/2 - 2e, 

(A7) 

where the minimal Abelian symmetry group is again Z%. 
• The constraints obtained from cases A 2 , A 3 , B\, and B 2 cannot be satisfied for any N. 



[1] G. C. Branco, L. Lavoura and F. Mota, Phys. Rev. D 39, 3443 (1989). 

[2] H. Fritzsch, Z. -z. Xing and Y. -L. Zhou, Phys. Lett. B 697, 357 (2011) [arXiv:1101.4272 [hep-ph]]. 
[3] G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher and J. P. Silva, arXiv:1106.0034 
[hep-ph] . 

[4] J. F. Gunion, R. Vega and J. Wudka, Phys. Rev. D 42, 1673 (1990). 

[5] W. Grimus, A. S. Joshipura, L. Lavoura and M. Tanimoto, Eur. Phys. J. C 36, 227 (2004) 
[hep-ph/0405016]. 

[6] J. L. Diaz-Cruz, J. Hernandez-Sanchez, S. Moretti and A. Rosado, Phys. Rev. D 77, 035007 
(2008) [arXiv:0710.4169 [hep-ph]]. 

[7] M. Frigerio, S. Kaneko, E. Ma, M. Tanimoto, Phys. Rev. D71, 011901 (2005). [hep-ph/0409187]; 
P. H. Frampton, S. Matsuzaki, [arXiv:0806.4592 [hep-ph]]; P. H. Frampton, T. W. Kephart, 
S. Matsuzaki, Phys. Rev. D78, 073004 (2008). [arXiv:0807.4713 [hep-ph]]; P. Ciafaloni, M. Pi- 
cariello, E. Torrente-Lujan, A. Urbano, Phys. Rev. D79, 116010 (2009). [arXiv:0901.2236 [hep- 
ph]]; P. Ciafaloni, M. Picariello, E. Torrente-Lujan, A. Urbano, Phys. Rev. D81, 016004 (2010). 
[arXiv:0909.2553 [hep-ph]]; T. Fukuyama, H. Sugiyama and K. Tsumura, Phys. Rev. D 82, 036004 
(2010) [arXiv: 1005.5338 [hep-ph]]; T. Fukuyama, H. Sugiyama and K. Tsumura, Phys. Rev. D 
83, 056016 (2011) [arXiv:1012.4886 [hep-ph]]; A. Aranda, C. Bonilla, R. Ramos and A. D. Rojas, 
Phys. Rev. D 84, 016009 (2011) [arXiv:1105.6373 [hep-ph]]. 

[8] A. Aranda, C. Bonilla and A. D. Rojas, Phys. Rev. D 85, 036004 (2012) [arXiv:1110.1182 [hep- 
ph]]. 

[9] G. C. Branco, D. Emmanuel-Costa and C. Simoes, Phys. Lett. B 690, 62 (2010) [arXiv: 1001. 5065 
[hep-ph]]. 



10 



[10] A. Aranda, C. D. Carone and R. F. Lebed, Phys. Rev. D 62, 016009 (2000) [hep-ph/0002044]. 
[11] A. Zee, Phys. Lett. B 93, 389 (1980) [Erratum-ibid. B 95, 461 (1980)]. 
[12] K. S. Babu and V. S. Mathur, Phys. Rev. D 38, 3550 (1988). 

[13] M. Fukugita, T. Yanagida, "Physics of neutrinos and applications to astrophysics," Berlin, Ger- 
many: Springer (2003) 593 p. 

[14] J. Schechter, J. W. F. Valle, "Neutrino Masses in SU(2) x U(l) Theories," Phys. Rev. D22, 2227 
(1980). 

[15] I. P. Ivanov, V. Keus and E. Vdovin, arXiv:1112.1660 [math-ph]. 
[16] A. J. Buras, hep-ph/9806471. 

[17] A. J. Buras, S. Jager and J. Urban, Nucl. Phys. B 605, 600 (2001) [hep-ph/0102316]. 
[18] A. Dedes and A. Pilaftsis, Phys. Rev. D 67, 015012 (2003) [hep-ph/0209306]. 



